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Abstract 

We use lattice Monte Carlo simulations to study non-perturbatively the tension, 
i.e. the free energy per unit length, of an infinitely long vortex in the three-dimen- 
sional U(l)+Higgs theory. This theory is the low-energy effective theory of high- 
temperature scalar electrodynamics, the standard framework for cosmic string stud- 
ies. The vortex tension is measured as a function of the mass parameter at a large 
value of the Higgs self-coupling, where the transition between the phases is contin- 
uous. It is shown that the tension gives an order parameter that can distinguish 
between the two phases of the system. We argue that the vortex tension can de- 
scribe the physics of long strings without lattice artifacts, unlike vortex network 
percolation. 
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1 Introduction 



During the early stages of its evolution, the Universe went through a series of 
phase transitions. Topological defects, such as cosmic strings, domain walls and 
magnetic monopoles, may have been created as a result of these transitions [1]. 
It is possible that the topological defects have observable consequences even 
in the present Universe. 

It is customary to view the topological defects as classical solutions obtained 
from the high-temperature effective potential of the system. However, the 
validity of this mean-field approach is uncertain, especially in gauge theories. 
Indeed, there is no symmetry breaking related to the Higgs mechanism [2], and 
in some cases the phases are analytically connected to each other [3]. Clearly 
the mean-field approximation does not work near the transition in such cases. 
It is therefore necessary to investigate how the inclusion of the full effect of 
fluctuations changes the picture of defect formation. 

To study the properties of string-like defects in gauge theories, it is most conve- 
nient to concentrate on the Nielsen-Olesen vortices [4] of the U(l) gauge+Higgs 
theory. Although this is much simpler than realistic theories containing cosmic 
strings, it permits a first-principles field theoretic study. In Ref. [5] an effective 
three-dimensional (3d) theory was constructed to describe the thermodynam- 
ical behavior of the full four- dimensional (4d) theory at high temperatures 
near the phase transition. The effective theory, the 3d U(l)+Higgs theory, is 
equivalent to the Ginzburg-Landau theory of superconductivity [6]. It contains 
fewer degrees of freedom than the original theory and its ultraviolet behavior 
is also better understood [7], which makes it much more suitable for numerical 
simulations. 

In discussing phase transitions, the first question is to find a proper variable 
to signal one. Although, in a gauge theory, there is no local gauge-invariant 
order parameter, the first order transition in the U(l)+Higgs theory can be 
localized by finding a discontinuity in, say, |0| 2 [8,9]. This does not help if 
the transition is continuous. Then a non-local order parameter may exist. For 
U(l)+Higgs one such is the photon mass [9] which vanishes identically in the 
symmetric (Coulomb) phase and is non-zero in the broken (Higgs) phase. 

For the U(l)+Higgs theory one may suggest further order parameters related 
to the existence of vortices. In fact, several different effective models have been 
proposed (see, e.g., [6,10]). These cannot replace a first-principles numerical 
study. In Ref. [11] the total density of thermal vortex- loop excitations was 
studied as an order parameter. However, this suffers from the difficulty of 
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performing a continuum limit: small vortices also appear as lattice artifacts. 

In scalar theories, it has been proposed that one should instead look at the 
density of infinitely long vortices, i.e., percolation [12] (see also [13]). However, 
as is well-known, this definition could be subject to lattice artifacts (for a 
review, see [14]). The basic problem is that "percolation" is a geometrical 
property, not directly related to the free energy of the system, and possibly 
sensitive to how the discretization needed for simulations is made. 

In this paper, we suggest a better alternative in the context of gauge theories: 
the tension of an infinitely long vortex, i.e., its free energy per unit length. We 
show how one puts a vortex on a lattice and how one measures its free energy 
using Monte Carlo simulations. We show that it really is an order parameter, 
vanishing in the V — > oo limit in the same phase as the photon mass. The 
result is given in continuum MS units, which makes it possible to compare it 
with the mean-field value [15]. The vortex tension thus defined is a physical 
observable and has a well-defined continuum limit. We suggest that the vortex 
tension is a good first-principles tool for discussing the physics of long vortices, 
and might perhaps also be used as an input in some phenomenological non- 
equilibrium estimates in the regime where mean-field theory does not work. 

The structure of the paper is the following. The U(l)+Higgs theory is defined 
and its properties are reviewed in Sec. 2. In Sec. 3 we describe how the vortex 
tension can be defined and measured, and in Sec. 4 we relate it to the phase 
structure of the theory. The details of the simulations, as well as their results, 
are given in Sec. 5. Sec. 6 contains the conclusions. 



2 The U(l)+Higgs theory 

In Ref. [5], the heavy degrees of freedom were integrated out from the 4d high- 
temperature U(l)+Higgs theory (see also [9,16]). The resulting 3d theory has 
the action 



where = d{Aj — djAi. Everything is expressed as dimensionless quanti- 
ties, by scaling with appropriate powers of e§, and the mass parameter y is 
renormalized in the MS scheme with the scale fx — l. 




(1) 
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Figure 1. The phase diagram of the 3d U(l)+Higgs theory. The vortex tension 
measurements in this paper are carried out in the region of continuous transitions 
(x = 2) for -1.6 < y < 0.5. 

The phase diagram of the theory (see Fig. 1) has been studied with lattice 
simulations in Refs. [8,9]. On the lattice, there are two alternative formula- 
tions, the compact and the non-compact one, but they are expected to have the 
same continuum limit in 3d. It was shown in Ref. [17] that in the non-compact 
lattice formulation there are two phases, the Higgs and the Coulomb phase. 
There is no symmetry breaking and no local order parameter, but the mass 
of the photon acts as a non-local order parameter: it is zero in the Coulomb 
phase and non-zero in the Higgs phase. This differs from the compact lattice 
formulation at any finite lattice spacing, and from the SU(2)+Higgs model, in 
which cases the phases are analytically connected [3,18]. 

At the value x = 2 that we will use in this paper, the transition between 
the two phases has been observed to be continuous, but quite different from 
standard symmetry breaking transitions [9] . At small x, a first-order transition 
was found in accordance with perturbation theory [19]. 

The field equations derived from the action in Eq. (1) have Nielsen-Olesen 
vortices as solutions when y < [4]. The vortex tension can be calculated in 
the mean-field approximation by solving the field equations numerically. The 
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result can be written in the form 

T MF = ^ = -^S(V2^), (2) 
L x 

where the function £ , with the value S(l) = 1, has been calculated numerically 
in, e.g., Ref. [15]. From there we interpolate that 8(2) fa 1.32, which value we 
use for comparison with our lattice results. 

To study the tension T non-perturbatively in the quantum theory, we dis- 
cretize the space and define the system on a lattice with lattice spacing a. We 
use the non-compact formulation, i.e. the link field is a real number instead 
of a phase angle. The lattice action corresponding to the continuum theory in 
Eq. (1) is 



x,i<j > G x,i 

+ [•••] E + JrE [<nx)0(x)] 2 , (3) 

x PG x 

where ^(x) = «j(x) + a_,-(x + ?) - a^x + j) - ^(x), C/,(x) = exp^a^x)], 
x = (xi,X2,x 3 ) with 1 < Xi < Ni, 

Pg = (4) 
a 

and the coefficient of the quadratic term is [7] 



[■■■] = 



ft 



6 + 



y 3.1759115(1 + 2x) 



2ttA 



(-4 + 8a; - 8x 2 )(\ogQ(3 G + 0.09) - 1.1 + 4.6a; 
16tt 2 /?2 



(5) 



Note how remarkably simple and analytic the lattice-continuum relation is. 

To give a gauge-invariant definition for a vortex on a lattice, we define for each 
link [11,20] (We choose the opposite sign here!) 



^(x,x+i) = - [«i(x) + 7(x + i) - 7(x)] 7r , 



(6) 



where 7 = arg0 and [X] n = X + 27mx such that [X]„. G (— 7r, n}. For links in 
the opposite direction, Y( x ,x-?) = _ ^(x-i,x)- F° r eac h closed curve C of links 
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Figure 2. If two vortices pass through one cell, the vortex tracing algorithm must 
decide how to connect the vortices, and this leads to an ambiguity in the length 
distribution of vortices. 

on a lattice, we then define the winding number n c as 



Z7r lec 



(7) 



In each configuration, the winding number gives the number of vortices passing 
through the curve C. Note that we count vortices with multiple winding as 
separate vortices on top of each other. 

Recently, there have been many attempts to interpret the phase transition 
in this and in related models as a percolation transition of vortices [12,13]. 
The percolation point y p is defined as the value of y above which one can in 
a typical configuration find a vortex that extends through the whole lattice. 
However, there are ambiguities in this procedure if two or more vortices meet 
in one cell (see Fig. 2), and it is not obvious that y p should coincide with y c , 
which is the transition point signalled by the non-analytic behavior in the free 
energy. For example, in the Ising model, these two transitions agree at d — 2 
[21], but (3 C « 0.95(3 P at d > 2 [22]. 

The idea behind the percolation studies is that as the transition is approached 
from below, the fluctuations can create larger and larger vortices, and at the 
transition point, the typical size of a vortex loop diverges. In other words, 
the free energy per unit length T of a vortex decreases to zero. To avoid 
the problems of the percolation approach, we have chosen here to measure T 
directly. Since our definition of T is based on the properties of the gauge field, 
this approach cannot be used in globally symmetric theories. 
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3 Inserting a vortex and measuring its free energy 

Free energies cannot be measured by Monte Carlo simulations, but changes of 
free energies can. To measure the free energy of a single vortex, we define a 
continuous vortex number (=winding number) m, measure the change of the 
free energy when m — > m + dm, < m < 1, and integrate these changes of m 
from to 1. The measurements can as well be done for m varying from to 
any integer. 

On the lattice, we call a field periodic, if it has the same value at x^ — iVj + 1 
as at Xi — 1. If we require the fields and a to be periodic, i.e. use periodic 
boundary conditions, then for any planar cross section (say, in the (x, y)- 
plane), the total flux $ = Z) P ianc Q; i2 and the winding n to t defined in Eq. (7), 
vanish in every configuration. The expectation value of an observable 0[a, 0] 
is then given by 

(0)= J VaV<j)0[a,<p]e- s[a ' (t,] , (8) 

periodic 

with S[a,4>] given in Eq. (3). 

We now define a non-periodic constant field a such that &2(Ni + 1, 1, rr 3 ) = 2n, 
and <Sj(x) = otherwise. The magnetic flux through the system in the z- 
direction, calculated from the background a, is $ = ^boundary ^« = 2tt. Note 
that there is no analogue to this in the compact formulation, where a is only 
defined modulo 2n. 

Let us now replace a by a + ma (m G Z) in the integrand of Eq. (8): 

(0) m = J VaV(j)0[a + ma,<P]e- s ^ +m& ^ . (9) 

periodic 

Note that a change of variables a + ma — > a transforms this to the form of 
Eq. (8) with non-periodic boundary conditions. Clearly, all the configurations 
in this integral have the total flux $ = 2ixm, and from Eqs. (6), (7) we notice 
that the total winding is n tot = m, i.e., the number of vortices going through 
the lattice is m. It can be shown that any configuration in which $ = 27rm 
and n to t — m and the physical quantities 0*0, (j)*{x)Ui{x)ip{x + i) and are 
periodic, can be gauge transformed to one that is included in the integral in 
Eq. (9). Therefore, Eq. (9) describes a periodic system in which the flux and 
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the winding number have been fixecQ 

Since L^(x) = 1 always, and a^-(x) = everywhere except 5i2(iVi, 1,2:3) = 2n, 
we see from Eq. (3) that 

S[a + ma, <p] = S[a, <p] + (3 G ^(27rma 12 (A^ 1 , 1, x 3 ) + 2vr 2 m 2 ). (10) 

^3 



From this expression it appears as if we had effectively inserted m vortices 
along the line (iVi, 1,2:3), 1 < £3 < N 3 . However, in Eq. (9) we integrate over 
all periodic field configurations, and thus all observables are translationally 
invariant. In fact, the precise form of a is irrelevant as long as its values are 
multiples of 2tt: only the total flux $ affects the results. 

Eq. (9) means that the free energy of a system with m vortices can be written 



as 



F m = -ln J VaV(j)e- s ^ +m& ^ . (11) 



periodic 



Only when m G Z, are the action and the observables periodic, and boundary 
effects are avoided. This is the well-known flux quantization condition. Still, 
F m is mathematically well defined for any m G JR., and the tension of the 
vortex can be written as (cf. Eq. (2)) 

T = ^ = -LJ ^ = -L / dm( i^M )m . (12) 

aJM 3 aJM-i J am aJM 3 J am 




Substituting Eq. (10) to Eq. (12) yields 



± 

T = 2vr 2 /? 2 J dmW(m 


1 

= 2vr 2 /3 2 J dm 



'^a 12 {N 1 ,l,x 3 )) m + 2m 

X3 



7CN 3 



(13) 



It is the quantity W(m) that we calculate with Monte Carlo simulations. 



7 There are also periodic configurations in which 7^ 27rntot> but their energy is 
assumed to diverge logarithmically with the area. Thus we neglect them here. 
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It is useful to inspect what kind of finite size effects one can expect for T. In the 
Coulomb phase, the massless photon gives rise to power-like finite size effects. 
Their magnitude can be estimated by considering the pure gauge theory, since 
the flux is expected to be distributed homogeneously in the whole system. In 
the absence of scalar fields, the flux quantization condition m G Z does not 
apply, and the system is translationally invariant at all m e R. Thus 

(ai 2 (iVi, 1, x 3 )) m = - ma 12 {N u 1, x 3 ) = 2nm - l) , ( 14 ) 

which yields 

o 



In the Higgs phase, all the fields are massive, and the finite size effects are 
exponentially suppressed. 



4 The relation of vortex tension to thermodynamics 



The behaviour of the vortex tension T can be directly related to possible phase 
transitions in the system. This is simply because a vortex carries a magnetic 
flux, and a magnetic field contributes to the free energy. Let us discuss this in 
some more detail. 

To begin with, note that a system with winding m corresponds to the magnetic 
flux density B z = 2Trm(3Q/N 1 N 2 . Thus, the free energy F(B Z ) of a system 
through which goes a non-vanishing magnetic flux, equals F m in Eq. (11), 
with m = a 2 NiN 2 B z /2i{. For symmetry reasons, F(B Z ) = F{\B Z \). In the 
thermodynamical limit, B z can have any real values. Hence the derivative of 
F(B Z ) becomes meaningful, and we can calculate it from the vortex tension, 
if we assume a repulsive interaction between the vortices, as is observed at 
large x. A system with a finite number m of vortices becomes namely infinitely 
dilute as the volume increases, yielding F m = mFi, and thus (we denote 
f(B z ) = l / ~ 1 F(5 2 ), the free energy density at a fixed flux density B z ), 



df(B z 



d\B z 



1 F(B Z ) - F(0) a 2 iViiV 2 T 

v ,i^o — m — = ~^r {Fl ~ Fo) = m (16) 
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Figure 3. If T > 0, there is a cusp in the free energy density f(B z ) at B z = 0. Thus 
the function g(H z ), defined in Eq. (17), has a flat region at — H c < H z < H c . 

where V = a 3 N 1 N 2 N 3 is the volume of the system. 

Let us denote by g(H z ) the free energy density of a system in which the flux 
is allowed to fluctuate and an external field H z is applied: 

oo 

g(H z ) = --\n[ dB z e- v ^- H ^ — > mm(f(B z ) - H Z B Z ) . (17) 

V J V-^oo B z 

— CO 



In the thermodynamical limit, g{H z ) is thus given by the Legendre transform 
otf{B z ). 

If T > 0, the derivative of f(B z ) is discontinuous at B z = (see Fig. 3). Its 
Legendre transform g(H z ) will therefore have a flat region, i.e. g(H z ) = g(0) 
when \H Z \ < H c = T/2n. This is precisely the Meissner effect, since (B z ) = 
—dg/dH z — 0, if \H Z \ < H c . It also implies that 



d 2 g 



= , when T > 0, 
h z =o > , when T = 0. 



(18) 



Thus, the values of T are related to the analytic structure (in particular, phase 
transitions) of the free energy. 

At x = 2, we expect T to behave continuously. It is non-zero in the Higgs phase, 
but decreases with increasing y, approaching zero at some value y = y c . At 
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larger values of y, the tension stays zero. If this turns out to be true, Eq. (18) 
shows that the free energy density g is non-analytic at y — y c . 



5 Simulations and results 

The Monte Carlo simulations of the theory in Eq. (3) at zero and non-zero 
values of m have been performed at the parameter value x = 2, at (3 a = 4. 
For this value of fie, the correlation lengths of massive modes like the scalar 
or the photon in the Higgs phase are large in units of the lattice spacing a. 
Indeed, for the values of y considered here, —1.6 < y < 0.5, even the shortest 
correlation length, which is the scalar one, is always larger than about 3a. 
Thus finite a corrections to the measured quantities are expected to be small. 

The major part of the numerical work has been performed with the help of 
Hybrid Monte Carlo (HMC) update algorithms, quite similar to those used in 
lattice QCD simulations. HMC updates were implemented separately for the 
gauge fields and for the scalar fields. Within both HMC algorithms a leapfrog 
discretization based upon 8 time steps into the molecular dynamics time di- 
rection was chosen. Acceptance rates were appropriately tuned. We measure 
our statistics in units of sweeps, where each sweep consists of two completed 
trajectories, one in the gauge and one in the scalar fields. The use of HMC 
updating schemes is favored by the fact that the usual checkerboard decom- 
position of degrees of freedom can be avoided and, therefore, the adaptation 
of the computer code to a parallel architecture is simplified. Most of our simu- 
lations have been executed on parallel computers. Some simulations have also 
been repeated on work stations with the use of different random number gen- 
eration algorithms and therefore different random number sequences. These 
runs serve as a valuable cross-check, and no inconsistencies were detected. 

We also experimented with over-relaxation updates of the gauge fields alone. 
As we failed to achieve a significant improvement in the statistical quality of 
the simulation, we will not report on the details of that approach here. 

We simulate symmetric cubic lattices with volumes N 3 . Typical lattice sizes 
range from N = 16 to iV = 32, which currently represents our largest system. 
The expectation value W(m) in Eq. (13) is evaluated by means of Monte 
Carlo simulations on a discrete set of 21 values m 8 partitioning the interval 
m = [0, 1] into 20 equally spaced subintervals. The integral of Eq. (13) was then 
calculated by the trapezoidal rule and it was checked that systematic errors 
introduced by the discretization are much smaller than statistical errors. A 
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Figure 4. W(m) on a 32 J lattice at x 



typical Monte Carlo result for the quantity W (m) is displayed in Fig. 4 for the 
32 3 lattice at y — — 1. Statistical errors on each individual rrii measurement 
can be as large as 40%, as is expected for an energy difference calculation 
on such large systems. We employed between 10000 and 20000 Monte Carlo 
sweeps on each of the individual rrii measurements. Statistical errors for the 
vortex tension T were then calculated using jackknife error calculation and 
statistical error propagation. 



In Fig. 5 we display our finite-volume data for the vortex tension T as a 
function of y in the whole interval of y considered on a selected set of 16 3 , 
20 3 , 24 3 and 32 3 lattices. One identifies three different regions. For y < — 1 
(Higgs phase), the vortex tension T appears to be finite and finite size effects 
appear to be under control. Note that at y = —1, data from several lattice sizes 
cluster around a common mean value. For values of y larger than about y = 
(Coulomb phase), T is non- vanishing in a finite system. However, it shows 
large finite size effects, which appear to be rather independent of y. There is 
also a region of intermediate values of y in which finite size effects are present 
and presumably not of a simple form. In this paper we will concentrate on 
those values of y that are clearly located in the Higgs or the Coulomb phase. 
A study of the intermediate y region would require finite size scaling methods 
beyond the scope of the present paper. 
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Figure 5. Lattice results for T function of y on 16 3 , 20 3 , 24 3 and 32 3 lattices 
at x = 2 and (3q = 4. The solid straight line is the mean-field result. Actually, only 
its slope is unambiguously determined. 

For two values of y, y = and y = 0.5, we investigated in detail the possible 
existence of a power-like finite size behavior T oc N~ 2 , as is characteristic of the 
pure gauge theory (see Eq. (15)). The observation of such a finite size behavior 
at the given values of y within the Coulomb phase signals the homogeneity of 
the system under the response of an applied flux. It is a strong argument in 
favor of the existence of a massless mode or alternatively, a vanishing vortex 
tension T. In Fig. 6 we display our finite N data for the vortex tension T 
as a function of N~ 2 at y — and y = 0.5. The lattice sizes considered are 
N 3 = 16 3 , 18 3 , 20 3 , 22 3 , 24 3 , 26 3 , 28 3 and 32 3 at y = 0. The same lattices except 
iV 3 = 26 3 have been considered at y — 0.5. 



Both data sets have been independently fitted with the two-parameter form 

T(N) = Ro^M + T(N = oo). (19) 



Deviations of the parameter R from unity measure deviations from a behavior 
anticipated in pure gauge theory, and the parameter T(N = oo) corresponds 
to our infinite volume extrapolation of the vortex tension T. The fit to all the 



13 




Figure 6. Finite size scaling data for T at x = 2, f3c = 4, y = 0.0 (left) and 
y = 0.5 (right), as a function of iV~ 2 . The solid line corresponds to a straight 
line fit as described in the text. The triangle corresponds to the infinite volume 
extrapolation. 

data, as indicated by the solid straight lines in Fig. 6, gives 



#o = 1.06(7), T(N = oo) = -0.019(32) at y = 0, (20) 
#0 = 0.97(6), T(N = oo) = +0.015(41) at y = 0.5. (21) 

We notice that R agrees with unity within errorbars. It thus appears that 
the possible renormalization of the coefficient in front of the N~ 2 behavior 
induced by the presence of the scalar field is negligible. Both fits exhibit rea- 
sonably small Xdof values, namely Xdof = 1-5 at y = and xiof = 1-4 at 
y = 0.5. Nevertheless, it can be noted that there are some data points which 
scatter. Simulations in the Coulomb phase turn out to be statistically rather 
demanding. From Fig. 7 we observe that the Monte Carlo simulation is better 
behaved (all errorbars cross the central value) in the Higgs phase, which is 
natural since there is no massless mode there. 

One may alternatively fix the parameter value T(N = oo) to zero in Eq. (19), 
and attempt to describe the data with 

T(N) = R 2 ^ (22) 
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where now any deviation of the parameter 77 from the value 2 signals departure 
from homogeneity caused by the presence of a vortex. We find 



77 = 1.94(14) at y = 0, 
77 = 1.93(12) at j/ = 0.5, 



(23) 
(24) 



and comparable values of the parameter R as in the previous case. These 
values clearly agree with the theoretical predictions for a Coulomb phase. 

In the Higgs phase, at y — —1.0, the vortex tension T was determined on 16 3 , 
20 3 , 24 3 and 32 3 lattices, see Fig. 7. Within the statistical errors the data are 
constant. A fit of the data to a constant value has xlof = 0-3 and gives 



Comparing this value with the mean-field result Tmf = 2.07 from Eq. (2), 
we observe a 17% difference, i.e., the mean- field result is close to the mea- 
sured value. One source for the discrepancy is that our result was measured 
at a finite lattice spacing a, while the mean-field calculation was performed 
in the continuum. One should also note that actually the mean-field calcu- 
lation only predicts the slope dT/ dy (see Fig. 5). Namely, the choice of 
the renormalization scale /i is arbitrary at the mean-field level. If we chose 
some other scale, the value of y would change by a constant term, y(fi) = 
y(l) + (4 — 8a; + 8x 2 )/(167r 2 ) ln/x In fact, changing from fi = 1 to fi = 0.2 
would bring the results in perfect numerical agreement. To remove this ambi- 
guity, corrections to the mean-field result would have to be computed up to 
2-loop level. 

Finally, we point out that since the infinite volume extrapolation of the vortex 
free energy at the values y = 0.0 and y = 0.5 is fully compatible with a 
vanishing vortex tension, and at y = —1 with a non- vanishing value, the 
vortex tension indeed behaves as an order parameter. 



6 Conclusions 

We have shown by numerical simulations that the vortex tension T defined in 
Eq. (12) acts as an order parameter at large values of x, where the transition 
between the Higgs and the Coulomb phases is continuous. The point y c , at 
which T vanishes, coincides with the true transition point, since the free energy 
density g was shown to be non-analytic at y c . 



T(N = 00) = 2.51(9) at y 



= -1 



(25) 
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Figure 7. Finite size scaling data for T at y = —1, x = 2 and [3g = 4, as a function 
of N. There are no visible finite size effects. The horizontal solid line corresponds 
to the infinite volume extrapolation T = 2.51(9), and the dashed lines indicate the 
error interval. 

At y < —1, the measured values of T agree with the mean-field result in 
Eq. (2) within errorbars, when the ambiguity of the renormalization scale in 
the mean-field approximation is taken into account. At y > 0, the system was 
observed to behave very much like pure gauge theory, i.e. the change due to 
the presence of a scalar field was very small. Outside the vicinity of y c , mean- 
field approximation thus seems to describe vortices well in both phases of the 
theory even at x — 2, where fluctuations are large. However, we expect that at 
y ~ y c , there will be large deviations from the mean-field behavior. The study 
of this interesting region requires significantly more computational resources 
than used here. 

We performed the simulations at a finite value of the lattice spacing a. Al- 
though we believe that the vortex tension T has a well-defined continuum 
limit, we expect that there are 0(a) errors, in particular in the Higgs phase. 
However, as of now, we have not removed them by a continuum extrapolation 
a -> 0. 

In cosmology, the time evolution of the vortex network after a phase transition 
is typically treated at the mean-field level [23]. Our results give confidence 
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to that approximation, at least deep in the Higgs phase. However, to see if 
modifications are needed to the picture of defect formation in gauge theories, 
a more detailed study is needed near the transition temperature to find out 
the role of the fluctuations. For a complete picture of cosmic string formation, 
one would also have to take the non-equilibrium effects into account. 

Another interesting way to extend this first-principles numerical study would 
be to consider configurations with multiple vortices. In the context of super- 
conductors, it has been suggested and also numerically verified in simplified 
models that the interactions between vortices give rise to new phases. It is 
then theoretically interesting and computationally demanding to explore, how 
the non-analytic structure in the multiple vortex case is related to the phase 
diagram at a vanishing magnetic flux density, Fig. 1. 
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